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Abstract 

We investigate the Kosterlitz-Thouless transition for hexatic order on a free 

fluctuating membrane and derive both a Coulomb gas and a sine-Gordon 

Hamiltonian to describe it. In the former, both disclinations and Gaussian 

curvature contribute to the charge density. In the latter, there is a linear 

coupling between a scalar field and the Gaussian curvature. We derive renor- 

malization group recursion relations that predict a transition with decreasing 

bending rigidity k. Using the sine-Gordon model, we show that there is no 

KT transition on a deformable sphere. 
PACS numbers: 05.70.Jk, 87.22.bt, 68.10.-m 
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A fluid membrane can be modeled as an ideal fluctuating surface whose energy depends 
only on its geometry. If the membrane is free from external tension, the dominant term in 
this energy is the Helfrich-Canham bending energy with bending rigidity n. Such a 
membrane with a molecular length a is crumpled at length scales beyond the persistance 
length, = ae^'^'^/^^, at which the renormalized bending rigidity passes through zero 0. 
Membranes with hexatic tangent plane orientational order are stiffer than fluid membranes 
A flat membrane has quasi-long-range (QLR) hexatic order at low temperature and 
undergoes a Kosterlitz-Thouless (KT) disclination unbinding transition to a disordered state 
0]. A free hexatic membrane at low temperature exists in a "crinkled" state with QLR 
hexatic order and a nonzero long- wavelength bending rigidity 

Free hexatic membranes are predicted p[ to undergo a KT transition from the crinkled 
phase to the crumpled fluid phase. To our knowledge, however, there has been no systematic 
treatment of this transition with complete renormalization group recursion relations for the 
hexatic stiffness K and disclination fugacity y associated with the KT transition in addition 
to the bending rigidity k controlling membrane height fluctuations. In this paper, we present 
an analysis of the KT transition on a fluctuating membrane using both the Coulomb gas 
model obtained directly from the hexatic Hamiltonian and the sine-Gordon model dual to 
the Coulomb gas model. A more detailed presentation will appear in a separate publication. 

Both local Gaussian curvature and disclinations contribute to the charge density of the 
coulomb gas. The bare hexatic rigidity K plays the role of the square of one of the charges 
in the Coulomb gas. The other charge is determined by K and the Liouville action 0. As 
in electrodynamics, neither charge undergoes a renormalization in agreement with Ref. [0, 
which found that K did not renormalize. We flnd, however, that fluctuations in Gaussian 
curvature reduce the dielectric constant even in the absence of dislcinations and lead to a 
decrease with decreasing n of the effective hexatic rigidity that controls correlations of the 
hexatic order parameter. As a result, we flnd, in contrast to previous treatments that 
decreasing n leads to a KT transition from the crinkled to the crumpled state as depicted 
m Fig. m . We obt ain recursion relations for K, y, and k that yield a KT transition and 
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that allow us to calculate quantities such as the persistence length in the crumpled phase. 
We find that ~ ^kt^'^^'^^^^ where ^kt ~ aexp(7r/6Vr — Tc) is the KT correlation length 
when K is large and the persistance length is larger than the C,kt- 

The sine-Gordon model we derive applies to an open membrane and to closed membranes 
(vesicles) of arbitrary genus h. As on fiat surfaces, this model has a scalar field (p with a 
gradient and a cos0 energy. The principal coupling to geometry occurs via an interaction, 
similar to the dilaton coupling of string theory P|, proportional to (pR with an imaginary 
coefficient, where R is the scalar or intrinsic curvature. Using standard renormalization 
group techniques for the sine-Gordon model we obtain the same recursion relations as 
we obtain using the Coulomb gas. We further show that shape fiuctuations on a sphere, like 
gauge fiuctuations in an infinite two-dimensional superconductor, generate a "mass" for 
and depress the KT transition temperature to absolute zero. 

Membrane coordinates in three space are specified by a vector R('u) as a function of a 
two-dimensional parameter u = {u^, u^). The metric tensor gab = da^{u) ■ (9;,R(£t), its inverse 
g"-^ , and the curvature tensor Kab = N ■ 9a9fcR, where da = d/du"" and N is the unit normal 
to the surface, can be constructed from R. The energy of a free membrane is described by 
the Helfrich-Canham Hamiltonian 



where g = det gab and H = g°'^Kab is twice the mean or extrinsic curvature. A membrane 
with hexatic order has an order parameter if) = |?/;|e^*^. A membrane with generalized p-atic 
order has an order parameter ip = \ip\e'^P^ invariant under 6^6 + (27r/p). The elastic energy 



where Aa is the spin connection satisfying '-y°-^daAb = S, where 7"'' is the anitsymmetric 
tensor with magnitude g'^^"^, and S = det = R/2 is the Gaussian curvature. The total 
long wavelength Hamiltonian for a p-atic membrane is then Ti = Ti^ + Tig. In the presence 




(1) 



associated with spatial variations of 9 for system with p-atic order is 




(2) 
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of disclinations, which have strength q = k/p where k is an integer, daO becomes singular: 
daO = daO + Va- The singular part satisfies 'j'^^daVb = n, where n = 2Txg^^^'^ Y,a Qa^iu — Ua) 
is the disclination density arising from disclinations of strength qa at positions Ua- Thus 
—Ag{va — Aa) = '~iadh{n — S) , where = DaD" is the covariant Laplacian operator on 
a surface with metric tensor Qab- Using this relation, Tie can be expressed as the sum of a 
regular part depending only on 6 and a curved space generalization of the energy of a 2D 
Coulomb gas: He = 'Hrcg + He with 

Hc = IkJ d'u^p{-Ag)-'p, (3) 

where p = n — S is the total "charge density" arising from disclinations and Gaussian 
curvature. Then, using J V9e~'^^"'<^ = e~^^^ where PHl = — 5'L/(967r) and 5*^ is the 



Liouville action |Tl[ 

Sl = J d\^R{-Ag)-'R, (4) 

we can write the partition function as 

Z{k, K, y) = Tr, J VWDOy^ e-'^^'^^^^o (5) 
= Tr, / PRe-^^«-^^^-^^=?/^. (6) 



Here Tie and Tic are understood to depend on the total number of vortices and on 
the position and charge of each vortex, Trv is the trace over vortices with approptiate 
factors of A^!, and PR is the measure for fiuid membranes, which includes the Fadeev- 
Popov determinant and Liouville factor [O. If we include only the lowest strength vortices 



with qa = ±l/p, then Tr^ can be defined via 

Trv= E ^N+-N.,px ^ L , f['Du^]N+['^u~]N_, 
N+,N- JV + .JV_. J 

where N± is the number of disclinations with strength ±l/p, x = 2(1 — h) is the Euler char- 
acteristic of the surface, and [Vu'^]]\f^ = d'^^tldi'^t)]^^'^ / '^'^ is integration measure 



for the positions of the N± charges with sign ±. The Kroneker S imposes the constraint 
that the total dischnation strength equal the Euler characteristic x- open surface, 

the long-range Coulomb potential leads to = A^_, and we can use Try with x = 0. We 
can convert the Coulomb gas model [Eq. (|^)] to a sine-Gordon model by representing the 
factor e~^'^'', via a Hubbard-Stratonovich transformation, as 



where the factor e'^^-^ assures that e ^^'^ is one when p = 0. Then, setting 6 



(l/27r) J^'^ due integrating over [Du~^]n^[Du and shifting — > p0/(27r) 



N,0 



uj, we 



obtain Z = J VRT>(f)e ^, where 

— ly d'^Uy/g cos (p — iX J d^u(j)S, (7) 

where 7 = / {Ati'^K) and X = p/ (27r). This effective action is valid for a surface of arbitrary 
genus. The coupling between and geometry is via a (pS term analagous to the dilaton 
coupling of string theory 0. 

Eqs. (H), (H), and (|^ provide us with three equivalent versions of the partition function 
for p-atic order on a fluctuating membrane, all of which can be used to study the KT 
crinkled-to-crumpled phase transition. We focus first on the Coulomb gas model on a nearly 
flat surface. The Monge gauge with -u = x = {x,y) and R(x) = {'x,h{x.)) provides the 
most natural description of such a surface. In this gauge, the Gaussian curvature is 5 = 
— ^{d'^6ij — didj)dihdjh to lowest order in h. Interactions among charges in a Coulomb gas are 
screened by other charges. In our Coulomb gas, there are two kinds of charges: disclination 
charge and Gaussian curvature charge. We are interested in the charge density p = n — S . 
The dielectric constant associated with this charge density is 

e-\q) = l-{(3K/q')CM, (8) 

where 



is the charge density correlation function. This expression for is identical to the 
corresponding flat space expression except that replaces p in Cpp. The effective 

or renormalized stifness is related in the usual |]13[ way to the zero q limit of e~^: 
I3Kr = limq^Q i3Ke~^. Cpp{q) is the complete p — p correlation function. It depends on all 
interactions in Eq. (|^), including the Liouville term (31-iL- We are interested in situations in 
which the disclination density, or equivalently is small, and we consider flrst the case y = 
for which there is only Gaussian curvature charge. Then Cpp{q) = Cgg{q), and PHl + PTi-c 
is equal to PHc{S = 0) with K replaced by PK' = PK — (127r)^^ in agreement with Ref. [0. 
The correlation function Cgg{q) can be expressed in terms of a polarization bubble P{q) as 
Cgg(g) = P{q)[l + {pK'/q^)P{q)]~\ To lowest order in {pKy \ C%{q) = P{q) with 

P(a) _ 1 1 M^fc |q X k|^ _ 3 

Thus, in the absence of disclinations, the renormalized rigidity is PK = PK—{3/327r){K/Ky. 
It is straightforward to verify that the correlation g{x.) = (cos(0(x) — ^^(0) — Jq dSaA"-)) 
between a spin at the origin and a spin at x parallel transported to the origin decays with 
|x| with an exponent determined by K^i rather than K: g{x.) ~ \x.\^'^'^^^^, which reduces 
when y = to \x.\^'^'^^^ . 

Next we consider the order y"^ terms arising from disclinations. To lowest order in Pk, 
the disclination density correlation function Cnn{<l) has the same form as it has in flat space: 
C°„(g) = g247]-3y2 j^oo ^^^3-27r/3ii:^ where a is the short distance cutoff. The cross correlation 
function Cnsil) is also of order y"^ and equal to PKCgg{q)C'^n{q) / q^. An expansion in powers 
of PK/ Pk then yields 

{PKn)-' = {PK)-' + 4nV -(-) ■ (10) 

Ja Qj \Qi / 



This is identical to the equation for Kr on a rigid flat surface but with {PK)^^ = 
{PK)^^ + (3/327r)(/5fi;)~^ replacing {PK)^^. It yields the KT recursion relations for PK and 
y. To these we must add the equation governing the renormalization of n. To lowest order 
in the temperture P~^ and y"^, the equation governing k is identical to that derived in Ref. 
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0. To leading order in /3 ^ and ?/, we can replace K by K{1) in the equation for k. Our 
recursion relations to lowest order in (3~^ and y are 



dl 

dy ( T^jSK^ 



dl \ p 



,2 



y (12) 



dl Att\ A(3k)' ^ ' 

We have also derived these equations directly from the sine-Gordon model using the pro- 
cedure or Ref. They are what one would have expected. The first two equations are 
exactly the KT equations, but with (3K replaced by f3K. Similarly, the equation controling 
j3n is identical to that derived using a momentum shell renormalization procedure in Refs. 
10] but again with (3K replacing [3K. Eqs. ( pT]) to (|13]) have a fixed line corresponding to the 
crinkled phase at ?/ = 0, = (3K /A independent of p. At (3K = (3K* = 2/7r, the system 
becomes unstable with respect to the formation of disclinations, and there is a KT transition 
to a fluid phase with unbound disclinations. For [I3K)~^ > {(3K*)~^, (3K decreases with 
/, and the ratio j3K / f3K tends to zero so that Eq. ([T3|) for f3K tends to that for a free fluid 



membrane. An altenative route to the derivation of Eqs. ([TT|) and (|1^) is to integrate out 
height fluctuations to produce an effective theory for 9 (or cj) in the sine-Gordon model) on 
a coarse-grained flat surface. This effective theory is identical to the usual flat space theory 
with K replacing K. 

In the vicinity of the KT critical point, y = 0, [3K = (3K* , I3k = /3k* = jSK* /A, we can 
set f3K = (3K*{1 — x), and = (3k*{1 — z). To lowest nontrivial order in x, y, z, Eqs. ( pUf) 
to ( |13|) become 

(15) 

The flow lines for these equations are shown in Figs. |T] and 0. In the xy plane (Fig. |^), 
flows are identical to the KT flows. In the xz plane (Fig. |]a), there is some curvature in the 



flows toward the fixed line resulting from couplings to y. In addition, flows are away from 
the critical line x = 0. Neither of these effects was observed in previous treatments ignoring 
y @]. Fig. |I|b shows flows in the {l3K)~^ — {(3k)~^ rather than the {I3K)~^ — {I3k)~^ plane. 
This shows that the crinkled phase exists only in a region near the origin with both {(3K)~^ 
and {(3k)~^ small. Thus, decreasing the bending rigidity k will lead to disordering to the 
fluid phase. This effect is not apparent in the treatment of Ref. f^. 

The persistance length in the fluid phase can be obtained with the aid of Eqs. ([14D 
to (plsl) . (3K reaches when x = 1, and does not become negative. We may, therefore, 
use Eq. to determine x{l) for / < Iq where x{Iq) = 1. The KT correlation length is 
^KT = ae'° = aexp[-K/by/T — TJ. For / > /q, we set x{l) = 1. Thus z{l) satisfles Eq. (|T5|) for 
/ < Iq and dz/dl = 3/2 for / > Iq, and 



where Zq = 1 - {(3k/ (3k*) and F{1) = dle-^^^'^x{l) with ^{l) = (3/2) J^dl{l - x{l)). The 
persistence length is then = ae' where z{l*) = 1. For large k, I* > Iq, and = 
(2g-p'('o)^^^g47rK/3T_ Pqj, smaller k and for PK near the critical point, /* may be less than Iq 
and ^p can be less than ^kt- 

The sine-Gordon model of Eq. (0) applies to surfaces of arbitrary genus. It can be used 
to investigate the effect of nonzero Gaussian curvature on the KT transition. On a nearly 
spherical surface of genus zero, we can parametrize the surface by its radius vector as a 
function of solid angle Q: = Rq{1 + i]{Q))er, where is the unit radius vector. To 

lowest order in t], 



where V is the gradient on a spherical surface of unit radius. Note that there is a linear 
coupling between and rj. In the case of a nearly flat membrane, this coupling was pro- 
portional to 0(V^/?.)^ and thus of higher order in normal displacments of the surface. An 




(16) 




(17) 
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effective theory on the reference spherical surface can be obtained by integrating out the 
height fluctuations. The resulting action is 

£ = 1^1 dfi0[-V' + {K/k)]<P+ ^J^^^^ J dficos0. 

The important feature of this model is the "mass term", {K/ K)(fP' that supresses fluctuations 
in and keeps the system in the ordered phase of the sine-Gordon model or, alternatively, 
in the disordered phase of the Coulomb gas model. If K/n > 1, the transition is supressed 
altogether, much as the KT transition is supressed in an infinite 2D superconductor. If on 
the other hand K/k, <^ 1, there can be an effective transition. In particular, when k, — * oo, 
the KT transition on a rigid sphere |jl5[ is regained. 
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(a) 7r/2p2 7^/2^2 

FIG. 1. Renormalization flows obtained from Eqs. (0) and (jl|) in (a) the (x, z) or 

(/?K)-i)-plane and (b) in the ((/3/s:)~\ (/3K)-i)-plane. (b) shows that increasing either 

{^K)~^ or (/?«;) ~^ leads to melting of the crinkled phase. 





FIG. 2. Renormalization flows Eqs. ( p!4| ) and (|T5|). (a) shows three-dimensional flows in y, 
and z. (b) and (c) show, respectively, flows projected onto the (y, z) and (y, x) planes, which are 
similar to those in the (y, x)-plane of a rigid flat surface. 
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